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We review our new method, which might be the most direct and efficient way for approaching the continuum 
physics from Hamiltonian lattice gauge theory. It consists of solving the eigenvalue equation with a truncation 
scheme preserving the continuum limit. The efficiency has been confirmed by the observations of the scaling 
behaviors for the long wavelength vacuum wave functions and mass gaps in (2+l)-dimensional models and (1+1)- 
dimensional a model even at very low truncation orders. Most of these results show rapid convergence to the 
available Monte Carlo data, ensuring the reliability of our method. 
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1. INTRODUCTION 

The main purpose of our work is to approach 
the scaling region and extract physical results 
by analytic calculations. For a lattice calcula- 
tion, a basic requirement is that for weak enough 
coupling, the dimensionless quantities should sat- 
isfy the scaling law, predicted by renormaliza- 
tion group equation. For SU(N C ) gauge theo- 
ries in 3 dimensions, superrenormalizabilty and 
dimensional analysis tell us that the dimension- 
less masses aM should scale as 

aM M 
g z e z 

For (2+l)-dimensional compact U(l) and (3+1)- 
dimensional non-abelian gauge theories, aM 
should scale exponentially as 



aM — > exp(—b/g 2 ). 



(2) 



If the calculated M data converge to a stable 
value, we can get an estimate for the mass. 

There have been various analytic methods 
available in the literature (for a review see 0). 
The main difficulty of the conventional methods 
(e.g. strong coupling expansion) is that they con- 
verge very slowly and very higher order 1/g 2 cal- 



culations are required to extend the results to 
the intermediate coupling region. Unfortunately, 
high order calculations are difficult in practice. 

Recently, we proposed a new method [j2|-|| for 
Hamiltonian lattice gauge theory. This method 
consists of solving the eigenvalue equation with a 
suitable truncation scheme preserving the contin- 
uum limit. Even at low order truncation, clear 
scaling windows for the physical quantities in 
most cases have been established, and the results 
are in perfect agreement with the available Monte 
Carlo data. Here we review only the work on 
U(l) 3 , SU(2) 3 and 2 dimensional a model, while 
that for SU(3) has been summarized in ||. 

2. THE METHOD 

The Schrodinger equation = en|fi) on the 

Hamiltonian lattice for the ground state 

= exp[R(U)}\0) (3) 
and vacuum energy en can be reformulated as 

^2{[E h [Ei, R(U)]] + [Ei,R(U)][E h R(U)]} 



2a 



(4) 
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To solve this equation, let us write R(U) in or- 
der of graphs G, hi , i.e., R(U) = Y, A R n(U) = 
J2 n ,i C n,iGn,i(U)- Substituting it to (g), we have 
the Nth order truncated eigenvalue equation 



N 
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(5) 



By taking the coefficients of the graphs in 
this equation to zero, we obtain a set of non- 
linear algebra equations, from which C n ,i are de- 
termined. The similar method applies to the 
eigenvalue equation for the mass and its wave 
function ||. Therefore, solving lattice field the- 
ory is reduced to solving the algebra equations. 

The lowest order graph is quite simple: 
Ri(U) = C 1A {U P + h.c.) The first term in (|) 
doesn't generate new graphs, but the second term 
does, i.e. 

[Ei, G ni (U)] <E R n (U) + lower orders, 



[Ei,G ni (U)][Ei,G n2 (U)] 6 R ni+n . 2 



+ lower orders. 



(U) 



(G) 



Two questions arise: 

1) Should all the new graphs generated by the sec- 
ond term in (|^) be taken as independent graphs 
of order n\ + nil For abelian gauge theories, the 
answer is yes. For non-abelian gauge theories, 
because of the uni-modular conditions J|-||,[13), 
there is a mixing problem not only for the graphs 
of the same order, but also for graphs of different 
orders. The classification for independent graphs 
is particularly complecate for SU(3). 

2) For ni + «2 > N, should we keep the lower or- 
der graphs in [E h G ni (U)][Ei,G n2 (U)]7 To pre- 
serve the correct limit, at Nth order truncation, 
one should to DROP all these graphs. This is the 
essential feature of our truncation scheme, which 
differs sufficiently from the scheme in M. 



There have also been some other truncation 
schemes proposed in One of their major 

problems is the violation of the long wavelength 
structure or continuum limit of the equation, and 
consequently the violation of the scaling law ([!]) 
or (||) for the physical quantities. 

Let's see further why the equation (||) should 
be truncated in the way suggested at point 2). 
The continuum limit of a graph G n ^(U) is 



G n ,i(U) = eVL4 n ,* tr(? 2 ) 



+ a 2 B n ^ tr(VFf 



(7) 



with T the field strength tensor and V the co- 
variant derivative. It has been generally proven 
Q that in the continuum limit the second term 
of (|5|) term should behave as 



[Ei,G ni {U)\[Ei,G r . 



t (U)} oc e 2 a e 



Tr{VT t 



(8) 



To preserve this correct limit, when the equation 
(||) is truncated to the Nth order, all the graphs 
created by [E u R ni (U)][E u R n2 {U)\ forni+n 2 < 
N must be considered. On the other hand, all 
the graphs created by this term for n\ + ri2 > N 
should be dropped, even there are lower order 
graphs. Otherwise the partial sum of the lower 
order graphs would make this term behave in a 
considerably different (wrong) way. 

3. RESULTS 

Once the coefficients C n .i are obtained by solv- 
ing (||), we can use (||) and (0) to compute the 
parameters /jq and ^.2 in the vacuum wave func- 
tion for the long wavelength configurations U Jl(| 

\Q) = exp[-fj, / d D ^ 1 x trT 2 



^2 I d°- 1 x tr{VTf 



(9) 



The results for fi$ and [12 in 3 dimensional 
SU(2) gauge theory are shown in Fig. 1. Em- 
pressively, nice scaling behavior is obtained even 
at N — 3, and the data for A/g 2 > 4 are in good 
agreement with the Monte Carlo measurements 
IpToh . The order N = 4 data are also included in 
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Figure 1. Parameters in the vacuum wave func- 
tion of SU(2) 3 model. The dashed lines show the 
mean values for the Monte Carlo data. 
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Figure 2. The mass gap of the SU(2) 3 model from 
three different methods. The dashed lines show 
the continuum limit of the Monte Carlo data. 



this figure. Although there are no big differences 
between the results at N = 3 and N = 4, higher 
order calculation seems necessary to ensure that 
the results would finally occur at the correct val- 
ues. 

Figure 2 shows our results for the mass gap in 
SU(2) 3 at different truncation orders. For com- 
parison, the results from the truncation method 
of Llewellyn Smith and Watson (LS-W) g and 
14th order series expansion [jllj are also included. 
Again, even at N = 3, with satisfaction of the 
scaling law, our results are in best agreement with 
the data taken from the continuum limit of the 
Monte Carlo data @. 

As is mentioned above, for non-abelian gauge 
theories, the uni-modular conditions lead to the 
existence of different choices for independent 
graphs. In SU(2), because of trU^ = trU p , all 
the disconnected graphs can be transformed into 
the connected ones, used as an independent set 
of the graphs. Our results in Figs. 1 and 2 are 
from such a choice, while the comparison between 
different choices (connected, disconnected p) and 
inverse) has been made in . Of course, one is 
free to choose arbitrary set of independent graphs. 
A criterion for a good choice is that it is conver- 
gent more rapidly to the continuum limit than 
other ones at lower order truncation. This has 
been obviously demonstrated for SU(3) 



The most intriguing scaling law is the exponen- 
tial scaling (Q). Before investigating QCD in 3+1 
dimensions, we would like to test our method in 
a (2+l)-dimensional compact U(l) model, which 
has many properties of the realistic theory. Here 
there is no ambiguity induced by uni-modular 
conditions in SU(N C ) theories. Because the na- 
ture of the abelian group greatly simplifies the 
calculations, we can easily write a program at ar- 
bitary orders. 
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wavelength vacuum state of compact U(l) 3 . The 
dashed line is the expected scaling law. 
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Figure 4. Relevant quantity for the asymmetric 
mass Ma of compact U(l) 3 . The dashed line is 
the expected scaling law. 



The results for U(l) 3 are shown in Figs. 3 and 
4 respectively. At relatively low truncation orders 
(comparing to 16th order series expansion [pT|), 
a scaling window has been seen around 1 / g 2 = 1 , 
and such a window becomes wider with the order 
N. Most impressively, there is an obvious ten- 
dency to converge to the scaling curve (dash line). 
This implies that the obtained physical quantities 
occur at their correct values. Fitting the N=6 
data in the scaling region, we get |H| 

l — = 3.1120 x l(T 2 e:cp(2.54(3)/.g 2 ), 
ag 

{M A agf = 365(73)e:rp(-5.0(2)/ ff 2 ), (10) 

or fi M A = 0.59(5). 

The non-linear (l+l)-dimensional a model is 
another interesting application of our method. 
According to the theoretical expectation, 

1 , 2i o 2 , 

Ma<x—exp( j - §- . (11) 

9 9 8tt 

Calculations at N = 5,6,7,8 have been carried 
out. Preliminary results are quite encouraging, 
and they are in reasonable agreement with the 
Monte Carlo data |l6| . To reach the asymptotic 
scaling region, higher order calculation seems nec- 
essary. 



4. SUMMARY 

The reason for the success of our method is 
that the truncated eigenvalue equation preserves 
the continuum limit. The results for (2+1)- 
dimensional models and (l+l)-dimensional a 
models are presented to support the efficiency 
and reliability of our method. In conclusion, 
the eigenvalue equation with a proper truncation 
scheme may be the most direct and efficient way 
for extracting the continuum physics. 
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